Revealing local failed supernovae with neutrino telescopes 
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We study the detectability of neutrino bursts from nearby direct black hole-forming collapses 
(failed supernovae) at Megaton detectors. Due to their high energetics, these bursts could be 
identified - by the time coincidence of N > 2 or N > 3 events within a ~ f s time window - from 
as far as ~ 4 — 5 Mpc away. This distance encloses several supernova-rich galaxies, so that failed 
supernova bursts could be detected at a rate of up to one per decade, comparable to the expected 
rate of the more common, but less energetic, neutron star-forming collapses. Thus, the detection of 
a failed supernova within the lifetime of a Mt detector is realistic. It might give the first evidence 
of direct black hole formation, with important implications on the physics of this phenomenon. 



The gravitational collapse of a stellar core is one of the 
most extreme phenomena in our universe. There, matter 
is pushed to its limits of density, and most of the energy is 
emitted by a non-electromagnetic form of radiation, the 
neutrinos, rather than in the final explosion (supernova) 
that often follows the collapse. 

Neutrinos are true tracers of core collapse. Due to 
their long mean free path, they give a direct image of the 
outskirts of the collapsed core. Furthermore, they are the 
only emission - together with gravitational waves - that 
always accompanies a collapse! Indeed, it is predicted 
that 10-20% of collapses directly generate a black hole 
PQ, with a brief and strong phase of neutrino emission, 
and no explosion [2^0]. For these failed supernovae, the 
star simply disappears from the sky, leaving the neutrino 
burst as a unique messenger of the event. 

At present, the detection of supernova neutrino bursts 
is still limited by long waiting times, as current detectors 

- of O(10) kt mass - can only capture the 1-3 bursts per 
century in our galactic neighborhood [5H7]. Upcoming 
Mt scale detectors will start to overcome the time barrier: 
for the common, neutron star-forming collapses (which 
have an accompanying explosion), they have a volume 
of sensitivity of 1-2 Mpc radius 6 , where about ~ 1 
collapse per decade is predicted [BJ. By applying a 10- 
20% fraction, this translates into about 1-2 detections of 
failed supernova bursts per century, still discouraging for 
an experiment lifetime of a few decades. 

In fact, however, certain factors enhance the de- 
tectability of a neutrino burst from failed supernovae. 
First, the higher neutrino luminosity and average energy 
of failed supernovae corresponds to a larger distance of 
sensitivity, a distance that happens as will be seen here 

- to be just enough to bring within the range of observ- 
ability several major, supernova-rich galaxies located 3-4 
Mpc away. This fortunate circumstance can boost the 
expected detection rate significantly, similarly to what 
was discussed for the diffuse supernova neutrino flux [5]. 
Furthermore, the shorter duration of a failed supernova 
burst (0.5-1 s) makes it easier to identify: the time coin- 
cidence of two neutrino events within ~ 1 s or so might 
be sufficient for discrimination against background. 



The fact that detecting individual neutrino bursts from 
failed supernovae is realistic, with a Mt detector, implies 
the potential to reveal - possibly for the first time - the 
direct collapse of a star into a black hole, with several 
implications on the physics of this transition, such as 
the rate of accretion of matter on the collapsed core, the 
equation of state of nuclear matter, etc.. Here we elabo- 
rate on the idea of the enhanced detection rate of failed 
supernova bursts, and discuss its implications. 

Failed supernovae (or direct black hole forming col- 
lapses, DBHFCs) are predicted to originate from stars 
with mass above M m i n <~ 25 — 40 Mq (with Mq the mass 
of the Sun) [J [5], corresponding to 9-22% of all collapsing 
stars (see, e.g., [S])- Numerical simulations [31 l4l I5HT2"] 
indicate that their neutrino burst lasts ~ 1 s or less, and 
has up to L ~ 5 • 10 53 ergs luminosity, due to the rapid 
contraction of the newly formed protoneutron star pre- 
ceding the black hole formation. The produced electron 
neutrinos and antineutrinos, v e and v e , have especially 
high luminosity, Lq c ~ Lqs ~ 10 53 ergs, due to the high 
rate of electron and positron captures on nuclei. Their 
average energy can reach E 0e ~ E os ~ 20 — 24 MeV. 

Due to oscillations in the star, the v e flux in a de- 
tector is an admixture of the unoscillated flavor fluxes: 
F s — pFg + (l—p)F®, where x indicates the non-electron 
species, v x =v^ u^, v r , v T , and p is the v e survival proba- 
bility [13]. Following [T2l [14] , we consider p = 0- 0.68 
and give results for the energy-independent, limiting case 
of p — 0.68, unless otherwise specified. We take the neu- 
trino fluxes from fig. 5 of [H] [3T], for the Shen et al. 
equation of state of nuclear matter. This set of flux and 
oscillation parameters maximizes F s [13], and so it is ad- 
equate to estimate the maximum potential of detection 
of failed supernovae. 

For comparison, we also model bursts from neutron 
star-forming collapses (NSFCs). These last 10-20 s and 
have typical parameters L ~ 3 • 10 53 ergs, Log ~ Lq x ~ 
0.5 • 10 53 ergs, E os ~ 15 MeV, E Qx ~ 18 MeV. The 
spectra of the produced neutrinos in each flavor typically 
have the form of a power-law times an exponential [15] . 
We restrict to the case in which p has the same value for 
the two collapse types [T2J Q3] . 



2 



Let us consider the response of a 1 Mt water Cherenkov 
detector [T6HT8] to a neutrino burst. The dominant de- 
tection reaction is inverse beta decay, v e +p —> n + e + , 
which we model as in [19] . The expected positron spectra 
for the two collapse types are shown in fig. [T] The higher 
energetics of a failed supernova is evident in the figure. 
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FIG. 1: Positron energy spectra at a 1 Mt water Cherenkov 
detector from a neutron star-forming forming collapse (thin 
curves) and a black hole-forming forming collapse (failed su- 
pernova, thick curves), at distance D — 1 Mpc. Dashed 
curves: p = 0; solid: p — 0.68. The Shen et al. equation 
of state is used for the failed supernova |12) . Integrated num- 
bers of events are also given for p = 0.68 and realistic energy 
windows of detection (see text). 



An experiment looks for inverse beta decay events 
within fixed time and energy windows designed to max- 
imize the signal to background ratio |20| . Typical time 
windows could be At = 10 s and At = 1 s for neutron 
star-forming collapses and failed supernovae [32] . The en- 
ergy windows are limited by background at low energy; 
a threshold of about 16 MeV in positron energy seems 
realistic [20]. The windows could be defined as including 
at least 80% of the events above this threshold: we find 
the intervals E e = 16 - 33 MeV and E e = 16-56 MeV 
for NSFCs and DBHFCs respectively. A neutrino burst 
is identified ( "detected" ) if N > N m i n ^2 — 3 events are 
observed in the energy window with time separation less 
than At. Note that the number of events due to a failed 
supernova increases with increasing p (fig. [TJ, i.e., with 
larger survival of the more luminous original v e compo- 
nent. It can be as large as fi{D) ~ 64(1 Mpc/Z?) 2 , up 
to five times larger than that from a NSFC. Therefore, 2 
(3) events are expected from a failed supernova as far as 
D ~ 6 Mpc (D ~ 5 Mpc). 

Given the "true" number of events, /i, the probabil- 
ity of detection of a burst (i.e., N > N m i n positrons 
observed) in the detector is given by the Poisson distri- 



P(N min ,D)= V ^M e -MW 

Z ✓ rrj] 



(1) 



-N„ 



[7] . It is shown in fig. j2|i) as a function of the distance D 
for N min — 2, 3. The figure confirms the expectation of 
a larger range of sensitivity to failed supernovae, with a 
probability of detection as large as 0.8 for D — D s ~ 4 — 
4.5 Mpc, which can be thus considered a typical distance 
of sensitivity. The corresponding distance for a neutron 
star-forming collapse is D s ~ 2 — 2.5 Mpc. 




FIG. 2: a): The Poissonian probabilities to detect N > 2 
and N > 3 events at a 1 Mt water Cherenkov detector for 
neutron star-forming collapses (lower curves, red) and black 
hole-forming collapses (upper curves, blue). The results of fig. 
[ljwith p = 0.68 are used, b): the rates of the two collapse 
types (the lower curves refer to failed supernovae) within a 
radius D from Earth, with their uncertainties. They are taken 
from [6] with a fraction /gg = 0.22 of failed supernovae. 

We now come to the key point of this work: how the 
increased distance of sensitivity allows to probe a region 
of high core collapse rate. Fig. [2]d) gives the nearby rates 
of the two types of collapses, Rbh and Rns, within a dis- 
tance D, with their uncertainty. They are derived from 
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the collapse rate in [B] (which is obtained from a catalog 
of galaxies [21] with conversion factors between luminosi- 
ties and core collapse rates [22 ), under the assumption 
of a constant, distance-independent, ratio Jbh = 0.22 of 
failed supernovae [33]. These rates are higher than the 
cosmological average [BJ, and actual supernova observa- 
tions favor an even higher rate [7J . Therefore, our results 
based on fig. ^p) are conservative. 

Fig. [2Jd) clearly shows the rapid increase of the rates 
between 3 and 4 Mpc, due to the presence of several 
galaxies (mainly IC 342, NGC 2403, M 81, M 82, NGC 
4945 [BJ), in this interval of distance. This is well 
within the range of sensitivity for failed supernovae, but 
only marginally accessible for the less luminous NSFCs. 
Within the typical distance of sensitivity, D s , fig. [2Jd) 
gives a rate of ~ 0.04 — 0.10 yr -1 for failed supernovae, 
and ~ 0.07—0.14 yr -1 for neutron star-forming collapses. 
The two rates are comparable, showing that the increased 
distance of sensitivity for failed supernovae compensates 
in part for their rarity. 

One can calculate the expected rate of detections of 
bursts from DBHFCs within a distance D [7J: 

R d B et H (N mm ,D) = J2 &RBHaP(N mmi D t ) . (2) 

i, Di<D 

The sum is over bins of distance, with Di < D, and 
ARbh,i is the failed supernova rate in each bin, so that 
Si d <d ARbha — Rbh(D). An expression analogous 
to eq. Q holds for the rate of detections of NSFCs, Rffig. 
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FIG. 3: The expected rates of detections of neutrino bursts 
that originate within a radius D from Earth, as functions of 
D, for neutron star-forming and black hole-forming collapses 
(upper and lower shaded regions). All parameters are as in 
fig. [2] for the case N > 2. 

Fig. [3] gives R B e H an d R'ns as functions of the dis- 
tance. Naturally, for each supernova type the detection 
rate follows the collapse rate for D <C D s ; it then flattens 
for larger distances, reflecting the suppression due to the 
small detection probability (fig. [2^,). This flattening oc- 
curs around 4 Mpc for NSFCs, and at ~ 8-9 Mpc for DB- 



HFCs. Depending on the normalization of the collapse 
rate, the detection rates for the two collapse types reach 
Rffs ~ 0.05 - 0.13 y^ 1 and R% ~ 0.04 - 0.11 yr" 1 . 

Thus, failed supernovae have a chance to be detected 
within the lifetime of an experiment. Due to their contri- 
bution, the total rate of burst detections could be twice as 
large as previously estimated, with a maximum of about 
2 detections per decade. 

Typically, expected detection rates are considered 
promising if they exceed the corresponding background 
rates, so that an observed burst can be attributed to a su- 
pernova with substantial likelihood. Assuming that cor- 
related events can be identified and subtracted [50], the 
background is given by accidental coincidences of uncor- 
rected events within the energy and time windows. By 
rescaling the SuperKamiokande measurements [23J [24] to 
a Mt mass, we find the rates of uncorrelated events to be 
A = 1855 yr _1 (A = 680 yr -1 ) in the energy window for 
DBHFCs (NSFCs). 

The rate of coincidence of two (three) such uncorre- 
lated events in the time window is (for A At <C 1) 0J2 — 
X 2 At (oj 3 ~ A 3 At 2 ) 25]. For failed supernovae (At = 1 
s) we find cj 2 — 0.10 yr -1 and uj 3 ~ 6.4 x 10~ 6 yr -1 . 
The same quantities for the NSFCs time window are 
u} 2 — 0.15 yr -1 and u 3 ~ 3.1 x 10~ 5 yr -1 . For both 
collapse types, the background doublet rate is compa- 
rable to or only slightly higher than the burst rate, so 
two observed positron events might be sufficient to claim 
a supernova detection, depending on the details of the 
experimental setup, and three events should give prac- 
tically certain identification. For a neutron star-forming 
collapse, the identification will probably be confirmed by 
the observation of the supernova explosion at telescopes, 
unless obscuration is substantial. For failed supernovae, 
one would have to rely entirely on neutrinos, or, possibly, 
on the concident detection of gravitational waves, or on 
establishing the disappearance of the star [2B] . 

Let us discuss how our results vary with the param- 
eters. Rates depend on Jbh as R B e H ^ Ibh an d 
R'ns ^ (1 — Jbh), so rescaling is immediate. The depen- 
dences on p and on the detector's mass, M, are described 
in fig. [4J which gives R d B e H within a 10 Mpc radius as a 
function of M, for different values of p and for the central 
curves in fig. [2]b) . For comparison, the background rates 
are shown; they depend on the mass as w 2 oc A 2 oc M 2 
and L03 oc A 3 oc A/ 3 . Expectedly, Rffjj increases with M 
and with p, due to the increase of the number of events 
(fig. [1} and therefore of the distance of sensitivity. Be- 
yond ~ 1 Mt of mass u>2 > RjiHi so a ^ least three events 
will probably be needed to establish detection. For a 
5 Mt detector like the proposed TITAND [27], we get 
Rg* ~ 0.10 - 0.16 yr" 1 for N mm = 3. 

Results also depend on the equation of state (EoS) of 
nuclear matter. For the softer Lattimer and Swesty EoS, 
the neutrino output of a failed supernova is somewhat 
less luminous and energetic, typically with Eq s ~ 20 MeV 
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FIG. 4: Solid lines: rates of background, to (red online), and 
of detected failed supernova bursts, Rbh (blue online), for 
sources within 10 Mpc distance, as a function of the detector 
mass, for N > 2 and p = 0,0.2,0.4,0.68 (lower to upper 
curves). Note that lj £ R% e H for M <: 0.8 Mt. Dashed lines: 
the same results for N > 3, for which the background rate is 
negligibly small (horizontal line). 



and L 0e ~ 0.5 • 10 53 ergs dUUl]. This translates into a 
reduced distance of sensitivity and therefore a lower rate 
of detections. Using the fluxes in [12], varying over the 
oscillation parameters and the local supernova rate, with 
M=l Mt, f BH = 0.22 and N mm = 2, we find i?f e f ~ 
0.016 — 0.045 yr _1 within 10 Mpc radius. This is close 
t° Rdet /bh/(1 — Jbh), as expected if the neutrino fluxes 
were the same in the two collapse types. 

Summarizing, the detection of a failed supernova at 
Megaton class neutrino detectors might be a realistic 
possibility, with a rate of detections reaching about one 
per decade, depending on the parameters. This is com- 
parable to the rate of the more common neutron star- 
forming collapses, and is due to the larger distance of 
sensitivity to failed supernovae, that includes several ma- 
jor, supernova-rich, galaxies. The short, ~ 1 s duration 
of a failed supernova burst might allow its unambiguous 
identification already with the coincidence of two inverse 
beta decay events within this time interval. 

Even with low statistics, the detection of a neutrino 
burst from a direct black hole-forming collapse will have 
profound implications. It might be the first observation 
of a different branch or core collapse, confirming its exis- 
tence, and giving information on the local rate of failed 
supernovae. This could be especially interesting in con- 
nection with the observed rate of bright supernovae being 
lower than expected [55], thus allowing for a substantial 
fraction of failed supernovae. It could also give the excit- 
ing opportunity to witness the formation of a black hole 
in real time, marked by the sudden truncation of the 
neutrino burst |29j . Considering the strong dependence 
of failed supernova neutrino bursts on the equation of 
state, conclusions about it might also be possible, with a 
high rate of DBHFCs bursts favoring a stiffer EoS. 
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